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In this paper we study double vector meson production in 77 interactions at high energies and 
estimate, using the color dipole picture, the main observables which can be probed at the Interna¬ 
tional Linear Collider (ILC). The total 7(Qi) -|- uiQV} —>■ Fi + V2 cross-sections for Vi = p, cf>, J/ip 
and T are computed and the energy and virtuality dependencies are studied in detail. Our results 
demonstrate that the experimental analysis of this process is feasible at the ILC and it can be useful 
to constrain the QCD dynamics at high energies. 
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I. INTRODUCTION 

There is an increasing interest in the construction of a high energy electron-positron collider [l| . The primary goal 
of this new facility will be to carry out precision measurements of electroweak physics, including the Higgs boson 
properties. An important byproduct of this program will be the study of high energy photon-photon collisions 0 and 
the continuation, at energies one order of magnitude higher, of the measurements performed at CERN-LEP, almost 
fifteen years ago. Photon-photon collisions are a very clean laboratory for the theory of strong interactions - Quantum 
Chromodynamics (QCD) - where we can test details of the QCD dynamics at high energies, such as the evolution 
both in virtuality (Q^) and in energy {1/x) (for a review see, e.g. Ref. [3 ). It has motivated the development 
of a large number of phenomenological studies in the last two decades [4l-llq. In particular, several authors have 
discussed the possibility of use the scattering of two off-shell photons at high energy in e+ e~ colliders as a probe of 
the parton saturation effects in the QCD dynamics, which are predicted to be present in the high energy regime [II- 
Although the experimental results on several inclusive and diffractive observables measured in ep scattering at HERA 
and hadron - hadron collisions at RHIC and LHC suggest that these effects are already observed in the energy regime 
probed by current colliders, these observations still need further confirmation. 

The state-of-art framework to treat QCD at high energies is the Color Class Condensate (CCC) formalism [T^, 
which predicts the saturation of the growth of parton distributions, with the evolution with the energy being described 
by an infinite hierarchy of coupled equations for the correlators of Wilson lines - the Balitsky-JIMWLK hierarchy 
(for recent reviews see [13). In the mean field approximation, this set of equations can be approximated by the 
Balitsky-Kovchegov (BK) equation [l^ . As emphasized in Ref. Q , in general, the applications of the CCC formalism 
to scattering problems require an asymmetric frame, in which the projectile has a simple structure and the evolution 
occurs in the target wave function, as it is the case in deep inelastic scattering. Therefore the extension of the BK 
equation to the calculation of the 77 scattering cross section is not a trivial task. In Ref. Q we have discussed this 
generalization in order to use the solution of the BK equation as input of our calculations of the total 7 * 7 * cross 
sections and photon structure functions, which were compared with the LEP data. In particular, in Ref. (3 we 
have improved the treatment of the dipole - dipole cross section, which is the main ingredient of the description of 
the 77 interactions in the dipole picture. Differently from previous phenomenological studies, which disregarded the 
impact parameter dependence, we have proposed an educated guess for this dependence and demonstrated that the 
LEP data can be described in this approach. The high energy behavior of the observables predicted in Ref. [3 is 
largely different from those obtained in previous studies. This conclusion motivates us to review the analysis of other 
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observables which could be measured at the ILC. One promising observable is double vector meson production in 
77 collisions, which has attracted the attention of several theoretical groups in the last years, with the cross section 
being estimated in different theoretical frameworks [IQ, as, for instance, the solution of the BFKL equation and 
impact factors at leading and next-to-leading orders. In this paper we will estimate the total 7 (Qi) + 7(Q2) —>■ 
cross-sections for Vi = p, (j), J/ijj and T considering the improved treatment of the dipole - dipole cross section and 
the energy and virtuality dependencies of the total cross sections will be analyzed in detail. Our analysis is strongly 
motivated by the fact that our knowledge about vector meson wave functions has improved considerably over the 
last years with the progress of phenomenological studies of vector meson production at HERA. As a consequence 
the main ingredients of our calculations are constrained by LEP and HERA data and hence our predictions for the 
ILC energies have only one free parameter - the slope parameter Hy, - which determines the t - dependence of the 
cross sections. The magnitude of this parameter for different combinations of vector mesons is still an open issue that 
deserves more detailed studies. 

This paper is organized as follows. In the next Section we present a brief review of the formalism, discussing in 
more detail the vector meson wave functions and the dipole - dipole cross section, which are the main inputs of our 
calculations. In Section im we present our predictions for the production of different combinations of vector mesons. 
In particular, the dependencies of the cross sections on the energy and photon virtualities are analyzed in detail. 
Finally, in Section ITVl we present our summary. 


II. DOUBLE VECTOR MESON PRODUCTION 
A. The cross section 


Let us review the main formulas of vector meson production in the color dipole picture (for more details see, e.g. 
Ref. [l 3 |)- The relevant scattering process is 7*7* V1V2, where Vi stands for both light and heavy vector mesons. 

At high energies, this scattering can be seen as a succession in time of three factorizable subprocesses (See Fig. [T|): i) 
the photons fluctuate into quark-antiquark pairs (the dipoles), ii) these color dipoles interact and, iii) the pairs convert 
into the vector meson final states. Using as kinematic variables the 7*7* c.m.s. energy squared s = IF^ = (p + 9)^, 
where p and q are the photon momenta, the photon virtualities squared given by Qi = —q^ and Q2 = —p^, and t, the 
squared momentum transfer, the total cross section for double vector meson production is given by 


cr (77 —>• Vi V2) 


dt 


(icr(77 —>■ Vi V2) 
dt 


1 da{'-fj —>■ Vi V 2 ) 
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V1V2 


dt 


^=o 
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( 1 ) 


where we have approximated the t-dependence of the differential cross section by an exponential with By^ y^ being 
the slope parameter. The imaginary part of the amplitude at zero momentum transfer A(s, t = 0 ) reads as 


ZmA( 7 *j* V1V2) 



X 


J dz 2 d'^r 2 'i'l^fi{z 2 , r 2 , Qj) 'itn%{z 2 , r 2 )crdd{xi 2 ,ri,r 2 ), 


( 2 ) 


where 'L''' and 'I'^* are the light-cone wave functions of the photon and vector meson, respectively. The quark and 
antiquark helicities are labelled by h, h, n and fi and reference to the meson and photon helicities are implicitly 
understood. The variable ri defines the relative transverse separation of the pair (dipole) and zi (1 — zi) is the 
longitudinal momentum fraction of the quark (antiquark). Similar definitions are valid for r2 and Z2. The variable 
X12 will be defined later. The basic blocks are the photon wave function, 4 ''*', the meson wave function, and 

the dipole-dipole cross section, add- 


B. Wave functions 

In the dipole formalism, the light-cone wave functions ^(z, r) in the mixed representation (r,z) are obtained 
through two dimensional Fourier transform of the momentum space light-cone wave functions 'I'^ ^(z, k) [^ . This 
subject has been intensely discussed in several references (see e.g. Refs. [ 2 ll - I^ L In what follows we present, 
for completeness, some of the main formulas. The normalized light-cone wave functions for longitudinally (L) and 
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FIG. 1: Double vector meson production in 7 * 7 * interactions at high energies in the color dipole picture. 


transversely (T) polarized photons are given by 

fN, 




4 ^ K-h(^ef2z{l-z)Q^^^ 
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Kojer) 

27T 


(3) 

(4) 


where = z(l — z)Q^ + TOj. The quark mass m/ plays the role of a regulator when the photoproduction regime is 
reached. The electric charge of the quark of flavor / is given by e e/. 

One simple way to model the vector meson wave function is to assume, following Refs. [HI, that the vector 

meson is a quark-antiquark state and that the spin and polarization structure is the same as in the photon case. The 
transversely polarized vector meson wave function is then given by 


4'wi.a=±i(^> “ ± \/‘2Nc - 


1 


hh,x=±iy^~^ ^-•‘^z{l-z) 

and the longitudinally polarized wave function is given by 


{ie^'^'-[z5h,±5-h^^ - (1 - z)5h,^5j,^^]dr + rnf5h,±5-h^±} 4>T{r, z). 
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mj-Vi 

Mvz(l — z) 




(5) 


( 6 ) 


where = (l/r)9i. + and My is the meson mass. The overlaps between the photon and the vector meson wave 
functions read then 


= e/e- 
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Meson My/GeV 

fv 

mf/GeV 

Nt 

Rf/GeV-'^ 

Nl 

N 

1 

> 

T 

9.460 

0.236 

4.5 

0.67 

2.16 

0.47 

1.01 

J/V- 

3.097 

0.274 

1.4 

1.23 

6.5 

0.83 

3.0 


1.019 

0.076 

0.14 

4.75 

16.0 

1.41 

9.7 

P 

0.776 

0.156 

0.14 

4.47 

21.9 

1.79 

10.4 


TABLE I: Parameters of the “Gauss-LC” vector meson wave functions. 


where the effective charge e/ = 1/3, 2/3, 1/3, or l/-\/2, for T, J/V', $, or p mesons respectively. The assumption that 
the quantum numbers of the meson are saturated by the quark-antiquark pair and that the possible contributions of 
gluon or sea-quark states to the wave function may be neglected, allows the normalization of the vector meson wave 
functions to unity. The normalization conditions for the scalar parts of the wave functions are then 


iVc _ (h 

27r Jo z‘^{l - 
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( 10 ) 


Another important constraint on the vector meson wave functions is obtained from the decay width. It is commonly 
assumed that the decay width can be described in a factorized way: the perturbative matrix element gg —>■ 7 * —>■ 
factorizes out from the details of the wave function, which contributes only through its properties at the origin. The 
decay widths are then given by 
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( 11 ) 

( 12 ) 


r =0 


The coupling of the meson to the electromagnetic current, fv, is obtained from the measured electronic decay width 

by 


T'K—>.e+e~ 


3Mv 


(13) 


We need now to specify the scalar parts of the wave functions, (j)T,Lir,z). Dosch, Gousset, Kulzinger and Pirner 
(DGKP) [2l| made the assumption that the longitudinal momentum fraction z fluctuates independently of the trans¬ 
verse quark momentum fc, where k is the Fourier conjugate variable to the dipole vector r. In the DGKP model one 
chooses 5 = 0 in Eqs. (0), (0 and m- The DGKP model was further simplified by Kowalski and Teaney [ 2 ^ , who 
assumed that the z dependence of the wave function for the longitudinally polarized meson is given by the short- 
distance limit of z(l — z). For the transversely polarized meson they set z) oc [z(I — z)Y in order to suppress 
the contribution from the end-points (z —>■ 0,1). This leads to the “Gauss-LG” wave functions given by: 

(/t(d z) = iVT[z(l - z)Y exp(-r^/ 2 i?y), (14) 

</l(g z) = Nlz{ 1 - z) exp(-r^/2i?|). (15) 


The values of the constants Nt,l and Rtx in Eqs. (ED and ED, determined by requiring the correct normalization 
and by the condition fv = fv,T = fv,L, are given in Table HI It is important to emphasize that this model allows to 
describe the HERA data and the recent LHC data for the exclusive vector meson photoproduction in hadron - hadron 
collisions (see, e.g. Refs. [25l - [^ L 


C. The dipole-dipole scattering cross section 

At lowest order, the dipole - dipole interaction can be described by the two - gluon exchange between the dipoles, 
with the resulting cross section being energy independent (See, e.g. Ref. [1^). Taking into account the leading 
corrections associated to terms oc log(I/a:), as described by the BFKL equation, implies a power-law energy behaviour 
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for the cross section, which violates the unitarity at high energies. These unitarity corrections were addressed in Ref. 
[ 2 ^, considering the color dipole picture and independent multiple scatterings between the dipoles, and in Ref. [s^ 
considering the Color Glass Condensate formalism. In the eikonal approximation the dipole - dipole cross section can 
be expressed as follows 

cr‘^‘^(ri,r2,T) = 2 J (fbAf{ri,r2,b,Y) (16) 


where A/'(t'i, r 2 , b,Y) is the scattering amplitude for the two dipoles with transverse sizes ri and r 2 , relative impact 
parameter b and rapidity semration Y. In Ref. Q A/” was assumed to be given in terms of the solution of the BK 
equation obtained in Ref. disregarding the impact parameter dependence, and we have proposed a model for 
the b dependence, which limits the range of impact parameters which contribute to the cross section. The basic 
motivation for this model is associated to the fact that although the unitarity of the ^-matrix {Af < 1) is respected 
by the solution of the BK equation (obtained disregarding the b dependence), the associated dipole - dipole cross 
section can still rise indefinitely with the energy, even after the black disk limit {Af = 1) has been reached at central 
impact parameters, due to the non-locality of the evolution. Consequently, a more elaborated model for the impact 
parameter dependence should be considered in order to obtain more realistic predictions for the dipole - dipole cross 
section. In Ref. Q we assumed that only the range b < R, where R = Max(ri, r 2 ), contributes to the dipole - dipole 
cross section, i.e. we assumed that J\f is negligibly small when the dipoles have no overlap with each other (b > R). 
Therefore the dipole-dipole cross section can be expressed as follows Q: 

nR 

a‘^^{ri,r 2 ,Y) = 2N{r,Y) (fb = 2TrR‘^N{r,Y), (17) 

Jo 

where N{r,Y) is the solution of the BK equation obtained in Ref. [3l| disregarding the impact parameter dependence, 
which we denote rcBK hereafter. The explicit form of a'^‘^ reads 

a‘^'^{ri,r 2 ,Y) = 2 TTrlN{r 2 ,Y 2 )Q{ri - r 2 )+ 2 TrrlN{ri,Yi)e{r 2 - n) , (18) 

where Yi = ln(l/xi) and 


- w^+Qr 


(19) 


As demonstrated in Ref. Q , using this model we can describe the LEP data for the total 77 cross sections and photon 
structure functions. 

For comparison, in what follows we also will present the predictions obtained using the phenomenological model 
for the dipole-dipole cross section proposed in The inclusion of these predictions in our analysis, allows us to 
estimate the theoretical uncertainty present in ILC predictions, as well as to make comparisons with existing results 
in the literature. The dipole - dipole cross section proposed in Ref. 0 is the following 


<^i‘!b{ri,r2,Y) = N{ri,r2,Y) 


( 20 ) 


with CTg’^ = (2/3)cto, where ctq is a free parameter in the saturation model considered, fixed by fitting the DIS HERA 
data. In the above equation N{ri,r 2 ,Y) = N{reB,Y = hi{l/xab)), where 




j.2^2 

2 ' 2 ^ab 

rf + r^ 


Q1+Q2 + -I- 4mf 

W^ + Ql + Ql 


( 21 ) 


Moreover, as in Ref. 0, we also will present the predictions obtained using the phenomenological model for the 
forward dipole scattering N{r, Y) proposed in Ref. [s^ and updated in [s^, which was constructed so as to reproduce 
two limits of the LO BK equation analytically under control: the solution of the BFKL equation for small dipole sizes, 
r <C \IQs{x), and the Levin-Tuchin law for larger ones, r ^ 1/Qs{x). In the updated version of this parametrization 
[ 3 ^, the free parameters were obtained by fitting the new HI and ZEUS data. In this parametrization the forward 
dipole scattering amplitude is given by 


N{r,Y) 



2(7.-t 


ln(2/rQa) ^ 


In^ {brQs) 

5 


for rQs{x) < 2 , 
for rQs{x) > 2 , 


( 22 ) 


where a and b are determined by continuity conditions at rQs{x) = 2, 7 ^ = 0.6194, k = 9.9, A = 0.2545, Qg = 1.0 
GeV^, xo = 0.2131 x 10“^ and A/q = 0.7. Hereafter, we shall call the model above HM-S. The first line from Eq. 
(E21) describes the linear regime whereas the second one includes saturation effects. One of the main motivations to 
use this model in our analysis is that it allows to estimate the magnitude of the saturation effects, by the comparison 
between the predictions of the full model with those obtained considering only the linear term. 
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FIG. 2: Energy dependence of the product Bv^v^ o'[7*(Qi)7* (Qi) V 1 V 2 ] assuming Vi = V 2 (Vi = p, </>, J/’I', T) and consid¬ 
ering Q? = Qi = 0 . 
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FIG. 3: Energy dependence of the product BviV 2 o'[ 7 *(Qi) 7 *(Q 2 ) —> E 1 V 2 ] assuming Vi / V 2 (V 1 V 2 = pJ/'^, , pT) 

and considering Qi = Q 2 ~ 0. 


III. RESULTS 

In what follows, as in Ref. Q, we will denote the predictions obtained using the dipole - dipole cross section given 
by Eq. (HOI) by model 1 and those using Eq. m as input by model 2. The parameters of our calculations are the 
same used in Ref. Q and this implies that our model gives a good description of the LEP data. As the value of 
the slope Bv^V 2 for th® different combinations of vector mesons in the final state is still poorly known, we will, in 
almost all cases, present our predictions for the product 14 V 2 ), which can be estimated without free 

parameters, since all parameters are constrained by the LEP and HERA data. 

In Fig. [2] we present our predictions for the energy dependence of the product BviV 2 o’ ( 7 *7* —t 1414) assuming 
14 = 14 (14 = P, 4)T) and considering = Qj = 0- In this case only the transverse photon polarizations 
contribute to the total cross sections. It is important to emphasize that the color dipole picture allows us to treat 
simultaneously double p production by real photons, which is a typical soft process, and double T production, 
which is the ideal laboratory to study the basic example of a hard process at high energies: the onium - onium 
scattering. Moreover, it allows us to study the transition between these two regimes, where we expect to see nonlinear 
(saturation) effects in the QCD dynamics. In our calculations we consider the two different models for the dipole - 
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FIG. 4: Dependence on the photon virtualities Qi = Ql = of fhe product BviV 2 {Q\)i*{O i) F 1 V 2 ] assuming V7 = V 2 
{Vi = p, 4>, J/^, T) for a fixed center-of-mass energy (W = 500 GeV). 


dipole cross section as well as the two models for the forward dipole scattering amplitude. We can observe that the 
main distinction is associated to the choice of the dipole - dipole cross sections. The predictions obtained using model 
2 are always smaller than those from model 1. Previous estimates of the double vector production have overestimated 
the magnitude of the total cross sections. This behavior was expected from our previous results for the total 7 * 7 * 
cross section Q. We also see that the difference between the predictions increases with the quark masses, going 
from a factor 4, in the pp case, to almost two orders of magnitude in the case of TT production. This suggests that 
the experimental analysis of double vector production at ILC can, in principle, constrain the model for the dipole 
- dipole interaction. Moreover, in the case of model 1, the IIM-S and rcBK predictions are almost identical for all 
combinations of vector mesons in the final state. In model 2 the IIM-S predictions are smaller than the rcBK ones for 
light vector meson production and larger for heavy vector meson production. Such behavior is directly associated to 
the distinct transition between small and large dipoles predicted by these two models for the forward dipole scattering 
amplitude (see Fig. 2 in Ref. 0 )- As expected, we find that our predictions are strongly dependent on the quark 
mass, with the cross sections being smaller for the production of heavier vector mesons. Similar conclusions are 
obtained in the analysis shown in Fig. [3j where we present our predictions for the energy dependence of the product 
BviV 2 <^[i*{QV) assuming Vi yf V 2 {V 1 V 2 = pJ/'k, J/'I'T, pT) and considering Qj = Ql = 0 . 

In Fig. Hlwe present our predictions for the dependence on the photon virtualities Qi=Q 2 = of the product 
BviV 2 {QV) ^ 1 ^ 2 ] for different combinations of vector mesons in the final state and fixed center-of-mass 

energy {W = 500 GeV). In this case we take into account the contributions of the transverse and longitudinal photon 
polarizations. For ^ 0 we have two hard scales present in the process: the mass of the quarks (vector mesons) 
and the photon virtualities. For the double light vector meson (pp, (jxf) production, the dominant scale is the photon 
virtuality. In this case our predictions strongly decrease with Q^. On the other hand, for the double T production, 
our predictions are almost - independent in the range considered, since the dominant scale that defines the size of 
the two interacting dipoles is the bottom quark mass. In contrast, for the double J/T production, the characteristic 
dipole sizes are determined at small by the charm quark mass and at medium by the photon virtualities. 
Consequently, we observe a mild dependence in the corresponding predictions. Moreover, we observe that the 
difference between model 1 and model 2 predictions increases at larger and for heavier vector mesons. In Fig. [S] 
we present our predictions for the production of two different vector mesons, which are similar to those observed in 
the production of identical vector mesons. Basically, the dependence is reduced for larger values of the sum of the 
masses of the vector mesons in the final state. 

In order to illustrate how the energy behavior depends on the masses of the final state mesons, on the photon 
virtualities Qi = Qi = ^.nd on the choice of the model for the dipole - dipole cross section, in Fig. [5] we present 
our predictions for the normalized cross sections. The different cross sections were all normalized to the unity at 
W = 100 GeV to better exhibit the different trends. For = 0 we observe a clear transition between the soft and 
hard regimes, with the growth with the energy being faster for heavier mesons in the final state. Moreover, we find 
that model 2 predicts a smaller slope than model 1. For = 20 GeV^, a similar behavior is observed, but in this 
case already for the pp production we see a steep rise of the cross section with the energy, which is directly associated 
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FIG. 5: Dependence on the photon virtualities Q\~ Q\ = of the product Bv^V 2 {Q 2 ) ^ 1 ^ 2 ] assuming V\ 7 ^ F 2 

(V 1 V 2 = J/^T,pT) for a fixed center-of-mass energy {W — 500 GeV). 
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FIG. 6 : Energy dependence of the normalized cross sections (see text) for different final states and different values of Q\ = 
gi = Q2. (a) g2 = 0 and (b) = 20 GeV^. 


to the presence of the hard scale . 

In certain cases, where the slope parameters are phenomenologically known, it is possible to make definite pre¬ 
dictions. In Fig. [7] we show the cross sections calculated with models 1 and 2 as a function of the energy W with 
the proper slope coefficients, taken from M- Bpp = 10 GeV Bpjj, = 5 GeV ^ and B^^ = 0.44 GeV Our 
predictions with model 1 are similar to those obtained in M , with small differences mainly associated to the different 
forward dipole scattering amplitude and to the treatment of the vector meson wave functions. In contrast, with 
model 2, we predict that at VF = 1 TeV the cross sections are ( 7(77 —>• pp) Ki 15 nb, < 7(77 —>■ pJ/'i') ~ 1.2 nb and 
t 7(77 — 7 > J/fk J/fk) Ri 0.25 nb, which are a factor r; 4 smaller than previous estimates in the literature obtained using 
the color dipole picture. 

Although saturation effects of QGD dynamics are expected to take over at higher energies, this change of dynamics 
can manifest itself with different strength for different observables. For this reason, when working with the dipole 
approach, which is naturally prepared to incorporate nonlinear corrections, it is always interesting to quantify the 
importance of the saturation effects. In Fig. [5] we show our results for the energy dependence of the 77 —> pp cross 
section, where we compare the full IIM-S model predictions with those obtained considering the linear regime of this 
model (first line in Eq. (I^^ b We see that the high ener gy behavior of the cross section is strongly modified by 
saturation effects. This conclusion was already obtained in [l4| and we see that it remains valid even after updating 
the dipole cross sections and model for the dipole-dipole interaction. 
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FIG. 7: Energy dependence of the 77 —^ hiV 2 cross section for different final states considering Qi = Q 2 = 0. 



W(GeV) 

FIG. 8 : Gomparison between the linear and full IIM-S predictions for the energy behavior of the 77 —>■ pp cross section. 


IV. SUMMARY 

The scattering of two off-shell photons at high energy in e’*' e~ colliders is an interesting process to look for parton 
saturation efffects. In these two-photon reactions, the photon virtualities can be made large enough to ensure the 
applicability of perturbative methods or can be varied in order to test the transition between the soft and hard regimes 
of the QCD dynamics. In recent years, a series of studies have discussed in detail the treatment of the total cross 
section and the exclusive production of different final states in 77 interactions considering very distinct theoretical 
approaches. One great motivation for these works is the possibility that in a near future 77 interactions may be 
investigated at the International Linear Collider (ILC). In particular, in Ref. Q we presented a detailed analysis of 
the 77 cross section at high energies using the color dipole picture and taking into account saturation effects, which are 
expected to be visible at high energies. In this paper we extended our approach to double vector meson production, 
improving the previous analysis in three important aspects: i) the theoretical treatment of the dipole - dipole cross 
section; ii) the forward scattering amplitude, considering the solution of the running coupling BK equation (which is 
the state-of-art of the CGC formalism); and iii) the treatment of the vector meson wave functions. Considering that 
all parameters of our approach have been fixed by fitting HERA and LEP data, our predictions for double vector 
meson production at ILC are parameter free, except for the only unknown parameter: the slope parameter BviVb, 
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which deserves a more detailed analysis. Our main conclusion is that the improvement of the theoretical framework 
for double vector meson production in 77 interactions resulted in a reduction of the previously estimated cross sections 
at ILC energies. However, our results indicate that the experimental analysis at ILC is feasible and may be useful to 
constrain the QCD dynamics at high energies. As a final remark we would like to say that understanding vector meson 
production in 77 collisions is very important not only for the phenomenology of future electron-positron colliders but 
also for exclusive double vector meson production in hadron - hadron collisions, which have been studied at the LHC 
and that could be further studied in future hadronic colliders. 
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